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Cover  photos  1  Very  finely  laminated  brownish  black  Grassy  Creek  Shale  with  very 
fine,  even,  and  closely  spaced  pyritic  laminae.  Large  black  feature  at  lower  left  is  a  pyrite 
nodule.  Sample  01KY01L1,  Christian  County,  Kentucky.  2  Finely  laminated  Grassy 
Creek  Shale  with  numerous  small  pyrite  nodules  along  bedding  planes.  Sample  02IL05C1, 
Effingham  County,  Illinois.  3  Poorly  laminated  Selmier  Shale.  Laminae  are  discontinu- 
ous and  widely  spaced.  Several  faint  silt  filled  burrows  cut  across  the  bedding  and 
truncate  laminations.  Sample  01KY08C1,  Christian  County,  Kentucky.  All  photos  are  at 
full  scale.  From  Cluff  et  al.  1981. 
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ABSTRACT 


Data  used  to  investigate  the  geochemistry  of  deposi- 
tional  environments  of  shales  of  the  New  Albany  Group 
(Devonian-Mississippian)  in  the  Illinois  Basin  were  the 
organic  carbon,  sulfur,  and  iron  contents  of  216  core 
samples  and  158  well  cuttings  samples,  reactive  iron 
determined  on  71  samples,  and  lithofacies  data  on  the 
core  samples.  The  positive  correlation  between  total  or 
pyritic  sulfur  and  organic  carbon  generally  found  for 
the  384  samples  indicates  that  the  dominant  process 
controlling  pyrite  formation  in  these  black  shales  was 
sulfide  production  by  bacterial  reduction  of  seawater 
sulfate  concomitant  with  bacterial  oxidation  of  organic 
matter. 

Shale  samples  with  a  low  organic  carbon  content 
(<2.5%)  generally  represented  dysaerobic  or  dysaero- 
bic/anaerobic  lithofacies  types  primarily  from  the  shal- 
lower western  and  west-central  parts  of  the  Illinois 
Basin.  These  samples  showed  a  good  correlation  be- 
tween sulfur  and  organic  carbon  (ratio  of  organic  carb- 
on to  sulfur  =  1.4  to  2.2)  and  between  degree  of 
pyritization  of  reactive  iron  or  total  iron  and  organic 
carbon.  Samples  with  greater  organic  carbon  content 
(that  is,  a  more  anaerobic  lithofacies  type)  were  gener- 
ally from  deeper  parts  of  the  Illinois  Basin  to  the  east 
and  south  and  had  organic  carbon  to  sulfur  ratios  aver- 


aging about  2.8.  Degree  of  pyritization  of  reactive  iron 
also  increased  and  total  iron  content  increased  slightly 
as  organic  carbon  content  increased.  Some  samples, 
mostly  from  the  Grassy  Creek  Formation,  had  rela- 
tively high  organic  carbon  content,  average  iron  con- 
tent, and  constantly  high  degree  of  pyritization  of  total 
iron;  pyrite  formation  was  probably  limited  by  the 
amount  of  reactive  iron  in  these  samples. 

Sixteen  organic  carbon-rich  samples  had  unusually 
high  iron  and  sulfur  contents.  Although  the  data  sug- 
gest pyrite  was  formed  in  H2S-bearing  bottom  waters 
in  five  samples  from  one  core  in  the  deepest  part  of  the 
Illinois  Basin  of  southeastern  Illinois,  anomalously  high 
amounts  of  pyrite  probably  formed  locally  during  later 
diagenesis  in  the  other  samples. 

Samples  from  the  top  of  the  Grassy  Creek  Shale 
from  the  deep  part  of  the  Illinois  Basin  (and  one  Hanni- 
bal Shale  sample)  had  relatively  high  organic  carbon 
content,  a  high  organic  carbon  to  sulfur  ratio  averaging 
4.8,  and  high  concentrations  of  vanadium  and  the  chal- 
cophilic  elements  Zn,  Se,  Sb,  and  Ag.  The  composition 
of  these  samples  reflects  slow  deposition  in  a  deepen- 
ing sea  when,  at  the  end  of  the  Devonian  or  beginning 
of  the  Mississippian,  the  southern  part  of  the  Illinois 
Basin  approached  sediment-starved  conditions. 


INTRODUCTION 


Investigations  of  the  relationship  between  organic 
carbon  and  pyritic  sulfur  in  marine  shales  have  led  to  a 
better  understanding  of  the  formation  of  pyrite  in  ma- 
rine sediments  and  of  the  depositional  and  diagenetic 
histories  of  many  shale  formations  (Leventhal  1983, 
Fisher  1986,  Leventhal  1987,  Donnelly  et  al.  1988,  and 
Boesen  and  Postma  1988).  Pyrite  forms  in  marine  sedi- 
ments from  the  reaction  of  iron  (usually  from  detrital 
minerals)  with  sulfide  that  is  formed  by  bacterial  reduc- 
tion of  sulfate  from  seawater.  Anaerobic  bacteria  use 
sedimentary  organic  matter  as  an  energy  source  and 
reducing  agent.  Sulfide  is  formed  only  under  anoxic 
conditions.  This  means  that  sufficient  organic  matter 
must  accumulate  in  the  sediment  for  some  of  it  to 
remain  after  all  the  available  oxygen  is  consumed  by 
aerobic  microorganisms,  which  convert  organic  matter 
to  carbon  dioxide.  Part  of  the  remaining  organic  matter 
undergoes  oxidation,  as  represented  by  the  simplified 
reaction 


2CH20  +  SCV 

organic  matter 


bacteria 


H2S  +  2HC03" 


When  sulfate  reduction  in  the  sediment  ceases,  remain- 
ing organic  matter  is  preserved  in  this  stage  of  diagene- 
sis. Organic  matter  may,  of  course,  be  decomposed 
further  via  fermentation  or  decarboxylation  reactions 
at  greater  depth  and  temperature. 


Under  normal  marine  conditions  of  oxygenated 
bottom  waters,  the  amount  and  bacterial  metaboliza- 
bility  of  organic  matter  appear  to  be  the  most  important 
factors  controlling  how  much  pyrite  may  form  (Berner 
1984).  Sweeney  (1972)  and  Goldhaber  and  Kaplan 
(1974,  fig.  10)  snowed,  for  modern  marine  sediments,  a 
correlation  between  organic  carbon  and  pyritic  sulfur, 
with  a  mean  ratio  of  organic  carbon  to  sulfur  (C/S)  of 
2.8  and  an  intercept  at  the  origin  (fig.  1).  The  line  of 
correlation  passes  through  the  origin  because,  if  no 
organic  matter  is  preserved,  there  is  no  net  reduction  of 
sulfate  to  sulfide  and  hence  no  formation  of  pyrite. 

In  euxinic  marine  environments,  such  as  the  Black 
Sea,  where  bottom  waters  contain  hydrogen  sulfide, 
the  metabolic  and  reduction  processes  of  the  bacteria 
operate  on  the  organic  matter  and  detrital  iron  sus- 
pended in  the  water  column.  Iron  sulfide  forms  in  the 
water  column  or  at  the  water-sediment  interface,  and 
an  appreciable  amount  of  pyrite  is  preserved  even  if 
little  or  no  organic  material  reaches  the  sediment  or  if 
the  organic  matter  is  largely  consumed  in  the  sulfate 
reduction  process  within  the  sediment.  Under  such 
conditions,  more  pyrite  is  formed  per  unit  of  buried 
carbon  (C/S  <  2.5),  and,  as  shown  in  figure  1,  a  regres- 
sion line  of  sulfur  on  carbon  gives  a  positive  intercept 
on  the  sulfur  axis  (Leventhal  1983). 

Raiswell  and  Berner  (1985),  using  data  from  recent 
Black  Sea  sediments  and  Jurassic  shales,  demonstrated 


Organic  carbon  (%) 

Figure  1  Schematic  plot  of  the  carbon  and  pyritic  sulfur 
relationships  in  modern  marine  sediments  and  in  samples 
from  a  euxinic  environment.  The  shaded  area  and  regression 
line  for  modern  marine  sediments  represent  the  bulk  of  the 
data  in  figure  10  of  Goldhaber  and  Kaplan  (1974).  The  shaded 
areas  and  regression  line  for  Black  Sea  sediments  represent 
data  of  Hirst  (1974)  as  interpreted  by  Leventhal  (1983,  fig.  1). 

that  the  limiting  factor  for  pyrite  formation  in  euxinic 
environments  is  the  amount  of  reactive  detrital  iron 
delivered  to  the  sediment.  Their  conclusion  is  demon- 
strated by  the  lack  of  a  relationship  between  organic 
carbon  and  the  iron  that  exists  as  pyrite.  Degree  of 
pyritization  of  iron  (DOP)  may  be  defined  (Berner 
1970)  as  (1)  the  ratio  of  pyritic  iron  to  total  iron  in  a 


sample  or  (2)  the  ratio  of  pyritic  iron  to  the  sum  of 
pyritic  iron  and  the  nonpyritic  iron  that  is  soluble  in 
hot,  concentrated  hydrochloric  acid.  The  second  defini- 
tion is  more  precise  because  nonpyritic  iron  that  is  not 
quickly  soluble  in  hot,  concentrated  hydrochloric  acid 
is  not  likely  ever  to  have  been  available  to  react  with 
hydrogen  sulfide  to  form  pyrite  (Berner  1970).  All  the 
reactive  iron  in  a  euxinic  environment  is  likely  to  be 
used  in  syngenetic  formation  of  iron  sulfide  in  the 
water  column  or  at  the  water-sediment  interface.  Con- 
sequently, the  degree  of  iron  pyritization  does  not  cor- 
relate with  organic  carbon  content  in  sediments  in 
euxinic  environments. 

Leventhal  (1987)  used  the  relationship  between 
organic  carbon  and  sulfide-sulfur  in  core  samples  of 
Devonian  shales  from  the  Appalachian  Basin  to  evalu- 
ate the  depositional  environment  of  the  shales.  He  con- 
cluded that  a  central  area  of  the  Appalachian  Basin  was 
euxinic,  as  implied  by  the  positive  intercept  on  the 
sulfur  axis  on  the  plot  of  carbon  versus  sulfur  content 
for  cores  from  that  area. 

This  report  investigates  the  relationships  between 
organic  carbon,  pyritic  sulfur,  and  the  degree  of  pyriti- 
zation of  iron  in  shales  of  the  New  Albany  Group 
(Middle  Devonian-Kinderhookian)  in  the  Illinois  Ba- 
sin. The  results,  together  with  lithofacies  data,  are  used 
to  characterize  environments  of  deposition  and  dia- 
genesis  for  these  shales.  This  report  also  examines  the 
relationship  between  these  environments  and  the  un- 
usual increase  in  concentration  of  trace  metals  in  the 
shales  at  the  top  of  the  New  Albany  Group  in  the  deep 
part  of  the  Illinois  Basin  in  Illinois.  Chemical  and  litho- 
facies data  are  taken  from  the  extensive  chemical  and 
geologic  studies  made  of  the  New  Albany  Group  in 
Illinois  for  the  U.S.  Department  of  Energy  Eastern  Gas 
Shales  Project  (Bergstrom  et  al.  1980).  In  addition,  re- 
active acid-soluble  iron  was  determined  in  71  New 
Albany  shale  samples. 


GEOLOGIC  SETTING  AND  STRATIGRAPHY 


The  Illinois  Basin  is  shown  in  figure  2,  and  the  forma- 
tions of  the  New  Albany  Group  in  Illinois  are  illus- 
trated in  cross  section  in  figure  3.  In  Indiana  and 
Kentucky,  and  recently  in  southern  Illinois  in  the  area 
of  the  Paducah  1°  x  2°  Quadrangle  (Nelson  1995),  the 
New  Albany  is  classified  as  a  formation.  In  the  remain- 
der of  Illinois,  it  is  classified  as  a  group.  The  group-level 
classification  is  retained  for  this  report  because  most  of 
the  data  are  from  the  parts  of  the  Illinois  Basin  outside 
the  Paducah  Quadrangle  in  Illinois. 

As  described  in  Cluff  et  al.  (1981),  the  basal  Blocher 
Shale  in  southeastern  Illinois  consists  of  brownish  black 
to  grayish  black,  finely  laminated,  dolomitic  and  cal- 
careous shales.  From  the  unit's  generally  high  organic 
carbon  content  and  the  finely  laminated  sediments 
with  rare  burrows  and  an  absence  of  megafossils,  Cluff 
(1980)  concluded  that  the  Blocher  Shale  was  deposited 
in  a  predominantly  anoxic  deep-water  environment. 


Overlying  the  Blocher  is  the  Selmier  Shale,  which 
consists  of  alternating  beds  of  thickly  or  poorly  lami- 
nated brownish  black,  grayish  black,  and  olive  gray  to 
olive  black  shales.  At  some  locations,  the  Selmier  com- 
prises indistinctly  bedded,  dark  greenish  gray  to  olive 
gray,  and  burrowed  shales  (Cluff  et  al.  1981). 

The  upper  part  of  the  Selmier  Shale  may  be  equiva- 
lent to  the  lowermost  part  of  the  Sweetland  Creek  Shale 
in  western  Illinois,  but  most  of  the  Sweetland  Creek  is 
younger  than  the  Selmier  (Cluff  et  al.  1981,  p.  28).  The 
Sweetland  Creek  consists  of  alternating  beds  of  dark 
greenish  gray,  greenish  gray,  and  grayish  green  shales 
that  are  indistinctly  bedded  and  moderately  biorur- 
bated.  Minor  amounts  of  olive  gray  to  olive  black  shales 
are  thickly  or  poorly  laminated  and  have  few  or  no 
burrows. 

The  Sweetland  Creek  Shale  of  western  and  central 
Illinois  interfingers  with  and  grades  laterally  into  the 
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Figure  2  Map  of  the  Illinois  Basin  with  locations  at  which  cores  and  well  cuttings  were  sampled  through  the  New  Albany 
Group.  Locations  are  reported  in  Bergstrom,  Shimp,  and  Cluff  (1980,  p.  12-14).  Paleogeographic  boundaries  are  from  Line- 
back  (1980)  and  Devera  and  Hasenmueller  (1991).  Locations  of  Eastern  Gas  Shales  Project  (EGSP)  cores  in  the  eastern  basin 
studied  by  Lechler  et  al.  (1979),  Shaffer  et  al.  (1981),  and  Hailer  et  al.  (1982)  are  denoted  by  D.  Locations  of  deep  basin  Grassy 
Creek  samples  and  of  a  Hannibal  Shale  sample  (from  07IL)  that  have  anomalously  high  trace  metal  contents  are  circled. 
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Figure  3     Generalized  northwest-southeast  stratigraphic  cross  section  through  the  New  Albany  Group  in  Illinois  (from  Cluff, 
Reinbold,  and  Lineback  1981). 


lower  part  of  the  Grassy  Creek  Shale  of  southeastern 
Illinois  (fig.  3).  Across  much  of  western  and  north-cen- 
tral Illinois,  the  Sweetland  Creek  Shale  is  not  consis- 
tently distinguishable  from  interfingering  and  over- 
lying olive  gray  to  olive  black  shales  of  the  Grassy 
Creek  that  have  thin  interbeds  of  grayish  olive  green  to 
grayish  green  shales.  But  in  southern  Illinois  and  west- 
ern Kentucky,  the  Grassy  Creek  Shale  is  a  sharply  de- 
fined unit  consisting  of  brownish  black  to  grayish 
black,  finely  laminated,  pyritic  shales. 

The  topmost  shale  units,  the  Saverton  and  the  Han- 
nibal, are  interbedded  greenish  gray  and  dark  greenish 
gray  shales.  The  units  have  some  olive  black  beds.  The 
greenish  gray  beds  are  nonlaminated  and  extensively 
burrowed.  The  Saverton  and  Hannibal  are  thickest  in 
western  and  central  Illinois  and  thin  to  the  southeast. 

On  the  eastern  side  of  the  basin,  in  south-central 
Indiana,  the  brownish  black,  calcareous  Blocher  Shale 
is  overlain  by  the  Selmier  Shale  (Hasenmueller  and 
Woodard  1981).  The  upper  part  of  the  New  Albany 
here  consists  of  three  members  that  are  not  always 


lithologically  distinct.  In  ascending  order,  these  are  the 
Morgan  Trail,  Camp  Run,  and  Clegg  Creek  Members. 
The  last  is  a  predominately  brownish  black,  finely  lami- 
nated, pyritic  shale.  These  member  names  are  also  ap- 
plied to  the  New  Albany  shale  strata  in  north-central 
Kentucky.  The  Morgan  Trail,  Camp  Run,  and  Clegg 
Creek  Members  are  time  equivalent  to  the  Grassy 
Creek  Shale  in  southeastern  Illinois  (Cluff  et  al.  1981). 

In  south-central  and  southeastern  Indiana,  several 
thin  beds  in  the  upper  Clegg  Creek  Member  (or  base  of 
the  Ellsworth  Member,  where  the  Ellsworth  Shale  can 
be  recognized)  are  formally  recognized.  In  ascending 
order,  these  are  (1)  the  Falling  Run  Bed:  a  0.2-foot-thick 
interval  of  phosphate  nodules;  (2)  the  Underwood  Bed: 
a  thin,  greenish  gray  shale  of  local  extent;  (3)  the  Henry- 
ville  Bed:  a  0.5-  to  3.5-foot-thick,  brownish  black  to 
black,  pyritic  shale  with  very  high  concentrations  of 
organic  carbon  and  trace  metals;  and  (4)  the  Jacobs 
Chapel  Bed:  a  thin,  greenish  gray,  glauconitic  shale 
(Hasenmueller  and  Woodard  1981,  Hasenmueller  and 
Leininger  1987). 


SAMPLES,  ANALYTICAL  METHODS,  AND  DATA 


For  the  Eastern  Gas  Shales  Project,  almost  400  samples 
from  the  New  Albany  Group  were  analyzed  for  organic 
carbon  and  sulfur  or  pyritic  sulfur,  as  well  as  for  all 
other  major  and  minor  elements  and  for  trace  elements 
(Frost  et  al.  1985).  Figure  2  shows  the  locations  of  the  15 
cores  from  which  216  samples  from  the  Blocher, 
Selmier,  Sweetland  Creek,  Grassy  Creek,  Saverton,  and 
Hannibal  Formations  were  taken  for  analysis.  In  addi- 
tion, 158  samples  (labeled  NAS-Onnn)  of  well  cuttings 
from  110  locations  (two  were  outcrops)  provided  cov- 
erage of  most  of  the  Illinois  Basin  within  Illinois  (fig.  2). 
For  the  Eastern  Gas  Shales  Project,  samples  were 
selected  from  the  cores  and  well  cuttings  as  follows. 
From  each  core,  a  3-  to  5-inch-long  segment  was  taken 
at  approximately  10-foot  intervals  through  the  New 
Albany  Group,  avoiding  pyrite  lenses  more  than  1  inch 
thick  or  any  particularly  bad  fracture  that  would  inter- 
fere with  sawing  out  the  segment  (Harvey  et  al.  1977). 
The  core  segment  was  split  vertically,  and  one-third  of 
it  was  crushed  and  ground  for  chemical  and  other 
analyses.  Sometimes  from  a  10-foot-long  core  section, 
an  additional  one  or  two  segments  were  sawed  and 
split  to  provide  samples  for  mineralogic  or  lithologic 
study;  these  samples  also  were  analyzed  for  chemical 
composition.  One  to  four  well  cuttings  samples  were 
analysed  from  any  one  well;  if  possible,  there  was  a  sam- 
ple from  every  100  feet  of  depth  and  at  least  one  sample 
from  each  formation.  Each  sample  was  picked  from 
cuttings  representing  a  5-  to  10-foot-depth  interval. 
Care  was  taken  to  avoid  any  obvious  caved  material. 

In  the  discussion  of  the  carbon-sulfur-iron  rela- 
tionships in  these  shales,  it  is  important  to  keep  in  mind 
that  the  samples  were  taken  at  wide  intervals  of  depth. 
Thus,  they  do  not  provide  a  continuous  or  inclusive 
record  of  carbon,  sulfur,  and  iron  contents  at  each  loca- 
tion through  the  formations  of  the  New  Albany  Group. 

Pyritic  Sulfur  Content 

Pyritic  sulfur  concentrations  were  determined  in  the 
216  core  samples,  but  not  enough  material  from  the  158 
well  cuttings  samples  was  available  to  do  analyses  for 
varieties  of  sulfur.  Results  for  the  core  samples  (Frost  et 
al.  1985)  showed  that  the  sulfur  is  nearly  all  pyritic. 
Pyritic  sulfur  averaged  94%  ±  6%  (one  standard  devia- 
tion) of  the  total  sulfur  in  all  core  samples.  The  percent- 
age was  highest  in  the  freshly  drilled  cores  and 
averaged  96%  ±  4%  in  the  samples  from  11  cores  drilled 
in  the  1970s,  compared  to  84%  ±  15%  in  34  samples 
from  4  cores  drilled  before  1970.  The  proportion  of 
sulfate  sulfur  was  highest  in  samples  from  core  01 IL, 
which  had  been  stored  in  air  since  1939.  These  results 
on  the  stored  cores  are  consistent  with  the  fact  that 
pyrite  undergoes  some  oxidation  to  sulfate  in  stored 
cores.  Therefore,  it  was  concluded  that  the  amount  of 
pyrite  originally  in  well  cuttings  samples  and  old  cores 
can  be  calculated  more  accurately  using  total  sulfur 
content  rather  than  the  pyritic  sulfur  content  that  is 


now  measured  in  a  sample.  In  this  report,  therefore, 
total  sulfur  concentrations  are  used  to  represent  pyritic 
sulfur. 

Organic  sulfur  is  ignored.  Organic  sulfur  values 
were  calculated  as  the  small  difference  between  total 
sulfur  and  the  sum  of  the  pyritic  and  sulfate  sulfur 
values  that  were  determined  for  180  of  the  core  sam- 
ples by  Frost  et  al.  (1985).  These  organic  sulfur  val- 
ues range  from  0%  to  0.6%  and  average  <0.1%;  they 
are,  however,  subject  to  large  absolute  analytical  error. 
Dickerson  (personal  communication  1979,  and  Dicker- 
son  1980)  calculated  an  average  organic  sulfur  content 
of  1.6  ±  0.3  weight  %  of  the  organic  carbon  content, 
based  on  elemental  analysis  of  the  benzene-methanol 
extracts  from  six  core  samples.  This  value  may  be  lower 
than  the  true  value  because  benzene-methanol  extrac- 
tion is  not  efficient  for  the  sulfur-containing  aromatic 
and  aliphatic  compounds  with  relatively  high  molecu- 
lar weights.  As  measured,  the  organic  sulfur  is  an 
insignificant  part  of  the  total  sulfur  and  insignificant 
with  respect  to  the  pyritic  sulfur  content  of  these  shale 
samples. 

Degree  of  Pyritization  of  Iron 

and  Determination  of  Reactive  Iron 

Elemental  concentrations  that  had  been  determined  in 
chemical  analyses  for  the  Eastern  Gas  Shales  Project 
were  used  to  calculate  the  proportion  of  the  total  iron 
residing  in  pyrite  (DOPt)  for  the  samples  for  this  re- 
port. DOPt  was  calculated  as  the  ratio  of  pyritic  iron 
(calculated  from  total  sulfur,  assuming  stoichiometric 
FeS2)  to  total  iron  determined  in  each  sample. 

As  previously  discussed,  the  proportion  of  the 
total  iron  residing  in  pyrite  may  not  reflect  the  fraction 
of  the  iron  that  could  have  formed  pyrite.  Some  iron  in 
sediments  might  occur  in  species  that  would  not  react 
with  hydrogen  sulfide  to  form  pyrite  (Berner  1970). 

To  test  whether  nonpyritic  iron  in  the  New  Albany 
shale  samples  studied  here  would  be  available  for  for- 
mation of  pyrite,  71  shale  samples  were  taken  for  deter- 
mination of  HCl-soluble  iron,  according  to  the  method 
of  Berner  (1970).  Replicate  0.1  g  portions  of  dried  sam- 
ples were  boiled  in  12  N  HC1  for  1  minute.  The  residue 
was  filtered  off  and  washed  with  distilled  water.  The 
HCl-leach  solution  was  diluted  to  50  ml  with  water. 
The  reactive  iron  concentration  was  determined  by 
atomic  absorption  spectrophotometry,  and  the  result- 
ing value  was  used  to  calculate  the  acid-soluble  iron 
content  of  the  shale  sample  (Appendix  1).  The  average 
mean  error  between  replicates  was  3.9%  relative. 

To  check  that  values  obtained  for  iron  as  deter- 
mined by  atomic  absorption  spectrophotometry  are 
comparable  to  those  obtained  by  X-ray  fluorescence 
spectroscopy,  the  same  71  shales  were  analyzed  for 
total  iron  by  atomic  absorption  spectroscopy,  using  an 
acid-digestion  procedure.  Resulting  values  (Appendix  1) 
agreed  with  the  total  iron  values  obtained  by  X-ray 


fluorescence  analysis  as  reported  in  Frost  et  al.  (1985), 
within  10%  relative,  and  with  random  difference. 

For  the  71  samples,  the  proportion  of  reactive  iron 
residing  in  pyrite  (DOPr)  was  calculated  as  the  ratio  of 
pyritic  iron  to  the  sum  of  the  pyritic  iron  and  hydro- 
chloric acid-soluble  iron  concentrations.  In  many  stud- 
ies, only  total  iron  in  a  shale  is  reported,  so  only  DOPT 
can  be  calculated.  To  evaluate  what  DOPj  can  reveal 
about  how  iron  controls  pyrite  formation  in  organic 
matter-rich  shales,  DOPj  values  of  New  Albany  shale 
samples  were  compared  to  the  DOPR  values. 

Lithofacies  Classification 

Lithofacies  data  from  previous  geologic  studies  of  the 
New  Albany  Group  for  the  Eastern  Gas  Shales  Project 
provide  independent  information  on  the  depositional 
history  of  these  shales.  Those  studies  used  the  stratified 
anoxic  basin  model  of  Rhoads  and  Morse  (1971)  and 
Byers  (1977),  which  correlates  the  upper  aerobic  zone 
(O2  >1.0  mL/L)  of  a  marine  basin  with  shelly  fauna  and 
a  rock  record  of  fossiliferous  limestones,  an  intermedi- 
ate dysaerobic  zone  (O2  =  0.1  to  1.0  mL/L)  with  biotur- 
bated  greenish  gray  mudstones,  and  a  bottom  an- 
aerobic zone  (O2  <  0.1  mL/L)  with  laminated  dark  gray 
to  black  shales. 


The  lithofacies  classification  for  the  shales  of  the 
New  Albany  Group  developed  by  Harvey  et  al.  (1977), 
Cluff  (1980),  and  Cluff  et  al.  (1980)  was  based  on  radi- 
ography of  shales  and  petrography  of  thin  sections  of 
samples  from  the  cores.  The  specimens  were  classified 
according  to  five  features  as  summarized  in  table  1. 
Because  the  lithofacies  classification  includes  color 
(which  is  largely  determined  by  organic  carbon  con- 
tent) and  pyrite  occurrence,  the  lithofacies  definitions 
do  not  describe  depositional  environments  in  a  way 
that  is  completely  independent  of  the  information  pro- 
vided by  the  organic  carbon  and  sulfur  concentration 
data  and  the  organic  carbon-sulfur  relationship  in 
these  shales.  However,  other  parameters  of  the  classifi- 
cation do  indicate  specific  new  aspects  of  a  sedimentary 
environment. 

At  times  during  the  formation  of  the  New  Albany 
Group,  aerobic  conditions  (>1.0  mL  02/L  water)  pre- 
vailed, and  limestones  and  dolomites  (for  example,  the 
Chouteau  Limestone,  fig.  3)  formed  in  shallow-water 
environments  (<50  m)  that  were  locally  high  in  energy 
(Cluff  et  al.  1980).  But  most  of  the  time,  oxygen  levels 
were  lower.  Then  the  shales  of  the  New  Albany  Group 
were  formed  in  depositional  environments  inferred 
from  the  four  lithofacies  types: 


Table  1     Lithofacies  characteristics  for  shales  of  the  New  Albany  Group  (after  Harvey 
et  al.  1 977  and  Cluff  et  al.  1 980). 


I 

II 

III 

IV 

Highly 

Indistinctly 

Thickly 

Finely 

bioturbated 

bedded 

laminated 

laminated 

mudstones 

shales 

shales 

shales 

Color 

gray  or  brownish  gray 

O 

greenish  gray 

X 

olive  gray 

X 

olive  black 

O 

X 

brownish  black 

O 

X 

Bedding  and  lamination 

massive 

X 

indistinctly  bedded 

X 

thickly  or  poorly 

laminated 

0 

X 

finely  laminated 

0 

X 

Bioturbation 

very  extensive  or  total 

X 

extensive 

X 

moderate 

X 

0 

sporadic 

X 

Pyrite  form 

fine  hair-like  features 

(burrows?) 

X 

X 

large  burrow  fillings 

X 

0 

nodules 

0 

X 

X 

laminae 

X 

X 

Spores 

absent 

rare 

common 

major 

X  =  major  or  common 


O  =  minor  or  rare 


Lithofacies  I  Bioturbated  mudstones:  dysaerobic 
(0.1  to  1.0  mL  O2/L  water),  moderately  deep  (50-150  m?), 
quiet  water  environments. 

Lithofacies  II  Indistinctly  bedded  shales:  mar- 
ginal dysaerobic  to  anaerobic  (<0.1  mL  02/L  water), 
moderate  to  deep  (150  m  or  greater?),  quiet  water  envi- 
ronments. 


Lithofacies  III  Thickly  laminated  shales:  anaero- 
bic, deep  (>150  m?),  quiet  water  environment,  broken 
only  by  occasional  currents. 

Lithofacies  IV  Finely  laminated  shales:  anaero- 
bic, deep,  very  quiet  water  environments. 

The  succession  from  lithofacies  I  to  IV  indicates 
qualitatively  increasingly  anaerobic  conditions. 


RESULTS  AND  DISCUSSION 


Organic  Carbon  to  Sulfur  Relationship 

The  organic  carbon  and  total  sulfur  contents  for  all  the 
New  Albany  shale  samples  are  plotted  in  figure  4, 
except  that  samples  from  the  Hardin  County  core  11IL 
that  would  plot  in  figure  4A  are  shown  in  figure  7.  Most 
of  the  samples  (fig.  4A)  show  a  general  positive  corre- 
lation between  organic  carbon  and  sulfur.  These  sam- 
ples plot  near  a  line  characterized  by  C/S  =  2.8,  which 
is  the  value  typical  for  sediments  from  modem  normal 
marine  environments  (Goldhaber  and  Kaplan  1974). 
The  few  remaining  samples  (fig.  4B)  have  a  much 
higher  or  lower  C/S  ratio.  Reasons  for  this  scatter  will 
be  discussed  below. 

Figure  5  shows  the  relationship  between  sulfur  and 
organic  carbon  content  for  core  samples  identified  ac- 
cording to  lithofacies  interpretations  made  by  Harvey 
et  al.  (1977),  Bergsrrom  and  Shimp  (1978a,  1978b,  1979), 
and  R.D.  Harvey  (personal  communication  1987).  The 
plot  and  the  values  in  table  2  indicate  that  organic  car- 
bon content  generally  increases  with  increasingly  anoxic 
conditions  (as  determined  from  lithofacies  interpreta- 
tion), although  there  is  considerable  overlap  of  points 
in  one  lithofacies  type  group  with  those  in  another. 

In  figure  5,  there  are  samples  with  carbon  and 
sulfur  values  that  plot  in  the  field  occupied  by  samples 
from  the  euxinic  Black  Sea  (Leventhal  1983,  fig.  1);  that 
is,  the  points  are  above  the  line  (shown)  drawn  by 
Leventhal  (1987)  to  separate  sediments  of  euxinic  envi- 
ronments (data  of  Hirst  1974)  from  sediments  of  nor- 
mal marine  environments  (data  of  Goldhaber  and 
Kaplan  1974).  Some  of  these  samples  that  plot  above 
the  line  are  low-carbon  content  samples  of  lithofacies  I 
and  II,  and  some  are  bioturbated  samples  of  lithofacies 
III.  For  example,  the  lithologic  description  of  the  sam- 
ple (S00174,  05IL01L3)  with  5.0%  organic  carbon  and 
3.5%  sulfur,  in  figures  4A  and  5B,  indicates  that  it  is 
from  a  thin,  bioturbated  bed  with  massive  laminations 
and  abundant  pyrite  nodules  3  to  7  mm  in  diameter 
distributed  in  burrows  (Bergstrom  and  Shimp  1978a). 
The  sample  in  figure  4A  (S00021,  01IL10L1)  with  1.3% 
organic  carbon  and  2.4%  sulfur  is  of  lithofacies  type  II 
(fig.  5B).  Moreover,  the  pyrite  in  the  sample  had  a 
positive  6  s,  which  indicates  sulfide  formation  from  a 
limited  reservoir  of  sulfate  (Anderson  et  al.  1987).  Sam- 
ples such  as  these  indicate  that  sediment  was  deposited 
in  at  least  partially  oxygenated  bottom  waters  and  that 
pyrite  formed  under  reducing  conditions  that  sub- 
sequently developed  in  the  sediments,  as  well  as  that 
some  pyrite  was  formed  in  a  later  stage  of  diagenesis. 


The  extrapolation  of  the  calculated  linear  regres- 
sion of  sulfur  on  carbon  for  various  groups  of  samples 
gives  a  positive  intercept  on  the  sulfur  axis  (table  2).  To 
test  whether  this  result  indicates  euxinic  conditions  in 
the  depositional  environment,  as  Leventhal  (1983)  re- 
ported, the  samples  from  core  06IL  were  evaluated. 

The  samples  from  core  06IL  (fig.  6A)  show  a  strong 
correlation  of  sulfur  with  organic  carbon,  except  for  one 
sample  that  has  unusually  high  sulfur  and  iron  con- 
tents. Excluding  that  sample,  the  10  samples  with  or- 
ganic carbon  contents  greater  than  2.5%  exhibit  a 
correlation  with  an  average  C/S  value  of  2.4.  Nineteen 
samples  with  relatively  low  organic  carbon  content 
(<2%)  have  an  average  C/S  value  of  1.4,  and  their 
preponderance  results  in  a  low  mean  C/S  value  of  1.7 
for  the  core  (table  2).  The  sulfur-carbon  covariance  in 
samples  from  06IL  is  well  fitted  by  a  parabola  drawn 
through  the  data  points  to  the  origin  (fig.  6A).  But, 
when  the  data  are  fitted  to  a  linear  regression,  the 
samples  with  low  organic  carbon  content  and  low  C/S 
ratios  contribute  significantly  to  causing  the  positive 
intercept  of  0.45  on  the  sulfur  axis  (table  2).  However, 
the  samples  with  organic  carbon  content  <2%  are 
mainly  of  lithofacies  type  II:  indistinctly  bedded,  bio- 
turbated sediments  with  pyrite-filled  burrows  (table  1). 
These  characteristics  indicate  that  the  pyrite  in  these 
samples  formed  later  in  diagenesis  after  the  sediment 
was  depleted  of  oxygen.  The  fact  that  pyrite  can  form 
at  several  stages  during  sedimentation  and  diagenesis 
complicates  the  interpretation  of  the  positive  y-axis 
intercept  of  the  sulfur  on  carbon  regression.  The  carbon- 
sulfur  pattern  shown  by  the  samples  from  core  06IL 
applies  to  the  other  groups  that  contain  many  samples 
with  low  organic  carbon  content:  core  04IL  (mean  C/S 
=  1.8  for  30  samples  [table  2]),  the  group  of  60  samples 
from  the  Saverton  and  Hannibal  Formations  (mean 
C/S  =  2.2  [table  2]),  the  lithofacies  I  and  II  groups,  and 
the  shallow  shelf  group  (table  2).  The  higher  level  of 
sulfur  relative  to  carbon  in  these  samples  may  be  due  to 
one  or  more  of  the  following  factors:  (1)  the  rate  of 
sedimentation  differs,  (2)  the  samples  reflect  a  source  of 
organic  matter  different  from  that  in  organic  matter- 
rich  samples  in  the  deep  basin,  (3)  aerobic  destruction 
of  organic  matter  is  greater  if  the  environment  cycled 
between  aerobic  and  dysaerobic  and  the  sediment  was 
reworked  by  the  burrowing  and  scavenging  activities 
of  benthic  organisms,  or  (4)  a  small  part  of  the  total  iron 
may  have  been  highly  reactive  so  that  iron  sulfide  ini- 
tially formed  easily  and  rapidly. 
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Figure  4  Variation  of  sulfur  with  organic  carbon  content  in  350  shale  samples  from  the  different  formations  of  the  New  Albany  Group 
in  the  Illinois  Basin.  The  straight  line  is  the  regression  line  (C/S  =  2.8)  for  modern  normal  marine  sediments.  The  dashed  line  is  that  drawn 
by  Leventhal  (1987)  as  representing  the  separation  between  sediments  from  euxinic  and  modern  normal  marine  environments.  (A)  Most 
of  the  samples  plot  in  the  field  about  the  regression  line  of  C/S  =  2.8.  Samples  from  core  11IL  in  Hardin  County,  Illinois,  that  would  plot 
in  this  field  are  shown  in  figure  7.  (B)  A  few  samples  tend  to  very  high  or  low  organic  carbon  to  sulfur  ratios.  Also  plotted,  from  Hailer 
et  al.  (1982)  and  N.R.  Shaffer  (personal  comm.  1988),  are  data  from  seven  samples  from  the  top  10  feet  of  the  New  Albany  (the  Clegg 
Creek)  in  cores  from  Marion,  Brown,  Jackson,  Clark,  and  Harrison  Counties  in  southeastern  Indiana  and  in  the  core  from  Sullivan  County 
in  southwestern  Indiana. 
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Figure  5  Variation  of  sulfur  with  organic  carbon  content  in  samples  classified  according  to 
lithofacies  type  from  cores  through  the  New  Albany  Group  in  the  Illinois  Basin.  The  straight  line 
is  the  C/S  =  2.8  regression  line,  and  the  dashed  line  is  Leventhal's  separation  of  euxinic  and 
modern  normal  marine  sediments. 
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Figure  5     continued. 
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Figure  6  Variation  with  organic  carbon  of  (A)  sulfur  and  (B)  iron  content  in  samples,  classified 
according  to  lithofacies  type,  of  New  Albany  Group  shale  from  the  Tazewell  County,  Illinois,  (06IL) 
core.  In  (A)  the  data  (excepting  one  datum  with  S  =  5.6%)  are  fitted  by  the  dashed  curve  whose 


.0:66 


function  is  y  =  0.72  x       for  organic  carbon  =  0.3%  to  5%. 
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Figure  7  Variation  of  sulfur  with  organic  carbon  content  in  25  samples  from  a  core  (11IL) 
from  Hardin  County,  Illinois.  The  samples  are  identified  (A)  by  formation  and  (B)  by 
lithofacies  type.  The  regression  equation  is  S  =  0.83  +  0.29  organic  carbon. 


Groups  of  samples  in  which  most  of  the  samples 
contain  more  than  2.5%  carbon  are  the  Grassy  Creek 
Shale,  the  Blocher  Shale,  the  samples  from  each  of  the 
shelf  margin,  basin  slope,  and  deep  basin  regions,  and 
the  samples  from  each  of  the  cores  from  those  regions. 
These  groups  have  higher  C/S  ratios  than  those  of  the 
low  organic-carbon  sample  groups  discussed  above 
(table  2),  with  mean  C/S  values  of  2.8  ±  0.4.  Excluding 
a  few  anomalous  samples  (the  samples  footnoted  in 
table  2  and  plotted  in  fig.  4B)  and  with  the  exception  of 
the  Blocher  Shale,  all  sample  groups  show  a  significant 


correlation  of  sulfur  to  carbon  concentrations.  Litho- 
facies and  carbon-sulfur  relationships  for  the  Illinois 
Basin  samples  are  consistent  with  a  oppositional  envi- 
ronment where  bottom  waters  are  dysaerobic  in  shal- 
low water  (where  organic-poor  sediments  are  depos- 
ited) and  become  increasingly  anoxic  in  quiescent 
deeper  water  (where  an  abundance  of  organic  matter  is 
produced  and  increasingly  preserved).  Iron-sulfur  and 
carbon-DOP  relationships  may  elucidate  in  more 
detail  the  paleoenvironmental  conditions  in  various 
regions. 
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In  table  2,  the  high  intercept  (0.83)  from  the  regres- 
sion of  sulfur  on  carbon  for  the  samples  from  core  11IL, 
Hardin  County,  Illinois  (fig.  7A),  requires  an  explana- 
tion different  from  that  advanced  for  the  other  cores 
discussed  above.  These  samples  have  a  wide  range  of 
organic  carbon  contents  (0.4%  to  6.5%),  a  change  to  a 
more  anoxic  lithofacies  type  with  increasing  organic 
carbon  content,  and  a  positive  correlation  of  sulfur  and 
DOPT  with  organic  carbon  content  (table  2).  However, 
the  samples  have  a  low  organic  carbon  to  sulfur  ratio 
(mean  1.6  ±  0.9)  over  the  entire  range  of  organic  carbon 
content,  thus  giving  a  high  intercept.  The  extra  pyrite  is 
particularly  evident  (fig.  7B)  in  a  few  samples  with 
<1.3%  organic  carbon  content  that  are  of  lithofacies  II 
classification.  The  low  C/S  ratio  may  be  due  not  to 
extra  diagenetic  pyrite  but  to  a  loss  of  organic  carbon. 
Relatively  high  temperatures  associated  with  sub- 
sequent igneous  intrusions  reported  from  the  Hicks 
Dome  area  where  this  core  was  taken  may  have  caused 
degradation  of  organic  material  and  loss  of  volatile 
components.  Evidence  for  thermal  effects  on  the  New 
Albany  shales  (depth  interval  0  to  261  ft)  in  core  11IL 
includes  negligible  release  of  gas  from  the  core  (Dick- 
erson  and  Chou  1980),  a  high  vitrinite  reflectance  gra- 
dient for  this  core  (Barrows  et  al.  1980),  a  lower  hydro- 
gen to  carbon  ratio  in  bitumens  in  samples  from  this  core 
than  was  observed  in  the  other  cores  (D.R.  Dickerson, 
personal  communication;  Frost  et  al.  1985,  p.  17)  and  a 
high  degree  of  crystallinity  of  illite  in  samples  from  this 
core  (Harvey  et  al.  1980). 


Iron,  DOPr  and  DOPT,  and  Organic 
Carbon  Relationships 

The  DOPr  values  for  52  samples  analyzed  for  HC1- 
soluble  iron  are  plotted  against  their  organic  carbon 
contents  in  figure  8.  Values  of  HCl-soluble  iron  were 
determined  and  DOPr  calculated  for  an  additional  19 
samples  from  a  group  of  metal-rich  samples  (fig.  10C). 
The  Blocher  Shale  samples,  which  are  high  in  carbon- 
ate, are  excluded  from  the  following  discussion  of  fig- 
ure 8.  The  strong  correlation  of  DOPr  with  organic 
carbon  at  low  levels  of  organic  carbon  (<3%)  is  associ- 
ated with  the  correlation  of  pyritic  sulfur  with  organic 
carbon  content  (fig.  4A)  and  indicates  that  formation  of 
pyrite  was  controlled  mainly  by  the  availability  of  reac- 
tive organic  carbon  for  bacterial  sulfate  reduction. 

At  higher  levels  of  organic  carbon  content  (>3%),  as 
organic  carbon  increases,  DOPR  continues  to  increase, 
but  with  great  variability,  although  the  correlation  be- 
tween the  two  parameters  remains  significant  at  the 
99%  confidence  level.  Some  of  the  variability  is  due  to 
the  wide  range  in  total  iron  content  of  these  shales, 
which  is  independent  of  organic  carbon  content.  This 
variation  is  shown,  for  example,  by  the  range  in  iron 
contents  of  several  samples  with  about  1%  organic 
carbon  content  in  figure  6B.  A  wide  variation  in  the 
detrital  iron  minerals  and  clay-adsorbed  iron  oxide 
contents  of  sediment  can  be  due  to  (1)  changes  in  the 
source  material,  (2)  the  grain  size  of  the  sediments, 
(3)  the  rate  of  production  of  organic  matter,  and  (4)  the 
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Figure  8     Variation  of  DOPr  with  organic  carbon  content  in  52  New  Albany  Group  shale  samples 
classified  by  formation.  The  Blocher  samples  were  excluded  from  the  data  fitted  by  the  curve. 
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rate  of  deposition  of  clastic  sediments  compared  with 
the  rate  of  deposition  of  organic  matter. 

Samples  from  core  06IL  also  show  a  slight,  but 
significant,  increase  in  iron  content  with  increasing 
carbon  content  (fig.  6B).  The  correlation  between  iron 
and  organic  carbon  content  is  not  significant  for  all 
groups  of  samples  as  grouped  by  core,  formation,  or 
lithofacies  type  (table  2).  Hence,  it  is  inferred  that  the 
main  source  of  iron  in  the  shale  is  detrital  material  or 
iron  adsorbed  on  the  clays  (Goldhaber  and  Kaplan 
1974).  Adsorption  of  colloidal  hydrous  ferric  oxide 
onto,  or  reduction  and  complexation  by,  organic  mate- 
rial provides  supplementary  means  of  accumulating  or 
holding  reactive  iron. 

The  data  points  for  the  71  New  Albany  shale  sam- 
ples that  were  analyzed  for  reactive  iron  are  scattered 
about  the  line  of  correlation  between  DOPr  and  DOPj 
(fig.  9).  This  scatter  of  data  points  suggests  that  some 
variable  amount  of  the  total  iron  was  never  available  to 
form  pyrite.  Therefore,  the  calculation  of  DOPT  does 
not  give  reliable  information  about  the  control  by  iron 
on  pyrite  formation.  At  higher  levels  of  DOPj  (fig.  9) 
and  organic  carbon  (fig.  8),  the  rate  of  increase  in  DOPR 
slows,  and  DOPr  no  longer  approaches  1,  which  may 
indicate  that  the  concentrated  HC1  treatment  overesti- 


mates the  amount  of  reactive  iron  (Middelburg  1991). 
The  DOPr  values  greater  than  0.85  (corresponding  to 
DOPj  values  greater  than  about  0.7)  probably  represent 
New  Albany  shale  that  lacks  significant  remaining 
reactive  iron.  That  is,  these  DOPr  values  represent  sam- 
ples that  are  iron-limited  with  respect  to  pyrite  forma- 
tion. This  representation  compares  with  an  average  for 
DOPr  of  0.84  in  two  systems  where  Raiswell  and 
Berner  (1985)  concluded  that  pyrite  formation  was 
iron-limited:  the  Lower  Jurassic  Jet  Rock  of  northeast 
England  and  the  Lower  Jurassic  Posidonia  Shales  of 
southern  Germany. 

Many  samples,  mostly  from  the  Grassy  Creek 
Shale,  have  relatively  high  organic  carbon  content  and 
DOPT  is  >0.7  (see  range  of  DOPT  in  table  2)  and  no 
longer  correlates  with  organic  carbon  content.  The  up- 
per limit  of  reactive  iron  needed  for  pyrite  formation 
apparently  was  reached  in  many  of  these  samples.  This 
iron  limitation  on  pyrite  formation  explains  why  some 
samples  with  high  organic  carbon  content  (>  -7%) 
show  no  increase  in  sulfur  content  (fig.  4). 

It  is  apparent  from  figure  8  that  pyrite  formation 
should  be  iron-limited  in  samples  with  very  high  or- 
ganic carbon  content.  Hence  it  is  surprising  that  a  cor- 
relation between  sulfur  and  organic  carbon  is  main- 
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Figure  9  Comparison  of  degree  of  pyritization,  DOPt,  calculated  as  the  ratio 
of  iron  in  pyrite  to  total  iron  in  the  sample  with  degree  of  pyritization,  DOPr,  as 
measured  by  the  ratio  of  iron  in  pyrite  to  the  sum  of  iron  in  pyrite  plus  HC1- 
soluble  iron  in  71  New  Albany  shale  samples. 
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Figure  10  Variation  of  (A)  iron,  (B)  DOPt,  and  (C)  DOPr  (18  of  the  samples)  with  organic  carbon 
content  in  22  samples  of  Grassy  Creek  Shale  that  have  unusually  high  trace  metals  content.  Samples 
are  from  the  upper  part  of  the  Grassy  Creek  from  15  locations  in  the  deep  basin.  Also  plotted  is 
Hannibal  Shale  sample  07IL01L2  (S00127)  that  has  the  same  type  of  high  trace  metals  content. 


tained  for  some  shale  samples  with  organic  carbon 
content  as  high  as  15.5%  (fig.  4A)  and  in  samples  with 
organic  carbon  contents  as  high  as  14%  from  several 
cores  from  south-central  Indiana  (Hailer  et  al.  1982). 
This  pattern  of  high  organic  carbon-sulfur  correlation 
may  be  an  artifact  arising  from  the  fact  that  in  the  set 
there  are  too  few  samples  high  in  organic  carbon  con- 
tent. To  the  extent  that  it  is  real,  one  possible  explana- 
tion for  the  correlation  is  that  pyrite  continues  to 
increase  with  increasing  organic  carbon  content  be- 
cause of  the  efficient  accumulation  of  reactive  iron.  For 
example,  calculations  on  the  35  samples  in  figure  8  that 


have  organic  carbon  contents  between  3%  and  12% 
show  correlations,  significant  at  the  95%  confidence 
level,  between  organic  carbon  content  and  total  iron 
(r  =  0.48)  and  between  organic  carbon  and  reactive  iron 
(that  is,  HCl-soluble  iron  plus  pyritic  iron  [r  =  0.66]). 
The  ratio  of  reactive  iron  to  total  iron  also  increases 
with  increasing  organic  carbon  content  of  the  samples. 
Also,  under  the  strongly  reducing  conditions  that  pro- 
duce organic  matter-rich  and  finely  laminated  shales  of 
lithofacies  type  IV,  iron-sulfide  reactions  continue  for  a 
long  time  period,  during  which  iron  that  is  less  reactive 
is  perhaps  also  employed  for  pyrite  formation. 
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Samples  with  Anomalous  Organic 
Carbon  to  Sulfur  Ratios 

Samples  with  low  organic  carbon  to  sulfur  ratios 

Random  high-pyritic  iron  content  is  shown  by  16  or- 
ganic carbon-rich  samples  (fig.  4B).  This  group  includes 
three  samples  with  the  lowest  organic  carbon  to  sulfur 
(C/S)  ratios  in  a  group  of  five  high-pyrite  samples  from 
Wayne  County  core  13IL.  These  five  samples  of  Grassy 
Creek  Shale  (S00412-S00416)  have  organic  carbon  con- 
centrations between  6%  and  9%,  very  high  sulfur  con- 
tents (>3.5%)/  very  high  iron  concentrations  (>4.5%), 
and  consistently  high  DOPt  values  (>0.65).  The  sulfur 
and  iron  contents  and  DOPt  values  do  not  correlate 
with  organic  carbon  content.  The  five  samples  are  from 
successive  10-foot  sections  (5,090  to  5,140  ft)  in  the  core; 
these  thinly  laminated  shales  show  no  evidence  of  bio- 
turbation.  The  high  pyrite  contents  in  samples  from  one 
50-foot  section  of  the  core  suggest  that  euxinic  condi- 
tions could  have  existed  locally  for  a  period  of  time 
during  deposition  of  the  Grassy  Creek.  Conversely,  the 
sulfur  isotope  composition  of  pyrite  in  these  samples 
does  not  show  sufficient  ^5  depletion  to  completely 
support  the  idea  that  euxinic  conditions  existed. 
Anderson  et  al.  (1987)  analyzed  the  pyrite  in  four  of  the 
samples  (S00412-S00414,  S00416)  and  found  that  6MS 
ranged  from  -17.6%o  to  -4%o.  These  values  are  more 
positive  than  the  values  expected  for  euxinic  condi- 
tions, that  is,  for  conditions  when  sulfide  forms  in  the 
water  column  and  where  access  to  ocean  water  sulfate 
is  unlimited.  For  example,  Beier  and  Hayes  (1989) 
found  d3^  values  of  -40%o  to  -30%o  for  shale  samples 
at  the  top  of  the  Clegg  Creek  Member  in  Indiana,  values 
that  indicated  to  them  syngenetic  pyrite  formed  within 
an  anoxic  water  column.  Therefore,  at  least  some  of  the 
pyrite  in  these  samples  from  core  13IL  formed  within 
the  sediments,  and  the  sediments  were  at  least  partially 
closed  to  seawater  sulfate. 


The  other  13  samples  high  in  pyrite  in  figure  4B  are 

(1)  five  samples  from  the  cores  06IL,  02KY,  01KY,  and 
11IL,  which  were  noted  in  table  2  and  one  sample 
(S00133  with  Fe  =  4.5%  and  S  =  3.5%)  in  core  07IL;  and 

(2)  seven  well  cuttings  samples  that  had  significantly 
higher  iron  (s~2o  above  the  mean  iron  level  for  sam- 
ples of  their  respective  formations)  and  sulfur  contents, 
including  one  Blocher  Shale  sample  (S00241  with  Fe  = 
5.3%  and  S  =  5.3%  on  a  carbonate-free  basis). 

These  13  samples  are  from  widely  scattered  loca- 
tions in  the  Illinois  Basin  and  from  various  formations. 
Specific  lithologic  information  is  not  available  on  the 
well  cuttings  samples.  Lithologic  characterization  of 
the  core  samples  (Bergstrom  and  Shimp  1978a,  1978b, 
1979;  Harvey,  personal  communication  1987)  indicates 
all  are  from  anoxic  environments  (lithofacies  III  or  IV) 
that  produced  shale  layers  with  even,  parallel,  and  thin 
laminations,  no  bioturbation,  and  abundant  small  py- 
rite nodules  distributed  either  randomly  or  along  bed- 
ding planes.  Thus,  the  high  pyrite  is  not  due  to 
anomalous  samples  that  contain  large  nodules.  A  high 
pyritic  iron  level  could  still  occur,  however,  from  sam- 
pling of  shale  with  pyrite  that  is  unevenly  distributed 
along  bedding  planes.  One  explanation  is  that  such 
erratic  pyrite  distribution  in  a  sample  reflects  the  local 
presence  in  later  diagenetic  stages  of  highly  reducing 
and  sulfide-rich  conditions,  which  could  result  from 
inhomogeneous  distribution  of  organic  material  or 
sulfate-reducing  bacteria  at  the  scale  of  microenviron- 
ments.  If  random  high  pyrite  content  does  indicate  a 
euxinic  environment  for  any  of  these  samples,  then 
such  an  environment  existed  only  within  a  narrow 
temporal  or  areal  distribution. 

Samples  with  high  organic  carbon  to  sulfur  ratios 
A  group  of  22  samples  from  15  locations  in  the  deep 
basin  (fig.  2)  that  are  characterized  by  higher  than  usual 
concentrations  of  trace  metals  (table  3)  generally  have 
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Table  3  Summary  of  elemental  concentrations  in  22  metal-rich  deep  basin  Grassy 
Creek  samples  and  in  metal-rich  Hannibal  Shale  sample  07IL01L2  (S00127)  com- 
pared with  concentrations  in  all  Grassy  Creek  and  undifferentiated  Grassy 
Creek/Sweetland  Creek  Shale  samples. 


185  Grassy  Creek 

22  deep  basin 

and  Grassy  Creek/ 

Grassy  Creek 

Sweetland  Creek1 
Geometric 

Geometric 

Range 

Hannibal  Shale 

mean 

(s) 

mean 

Sample  07IL01L2 

Organic  C  (%) 

7.21 

5.80-8.95 

6.1 

5.98 

Total  S  (%) 

1.58 

1.16-2.14 

2.3 

0.44 

Fe  (%) 

3.03 

2.58-3.56 

4.0 

3.41 

V  (ppm) 

426 

324-560 

210 

1,600 

Cr  (ppm) 

105 

85-127 

82 

180 

Co  (ppm) 

13 

11-16 

23 

30 

Ni  (ppm) 

115 

86-152 

84 

260 

Cu  (ppm) 

82 

66-103 

79 

60 

Zn  (ppm) 

403 

218-744 

150 

1,100 

As  (ppm) 

32 

23-44 

31 

6.8 

Se  (ppm) 

14 

10-20 

7  (38)2 

12 

Mo  (ppm) 

43 

20-92 

60 

56 

Ag  (ppm) 

0.9  (12)2 

0.4-1.9 

<0.2  (55) 

<0.8 

Sn  (ppm) 

6.7 

5.3-8.6 

6.6(180) 

17 

Sb  (ppm) 

8.0 

6.0-1 1 

3.9 

6.9 

Pb  (ppm) 

37(17) 

18-37 

37(155) 

42 

U  (ppm) 

22 

17-30 

26 

16 

From  Frost  et  al.  (1 985,  p.  23).  Samples  from  core  1 1 1L  not  included 
No.  of  samples  with  real  values 


high  organic  carbon  to  sulfur  ratios.  (The  22  samples 
are  identified  and  the  elemental  concentrations  are 
listed  in  Appendix  2.)  All  of  these  samples  are  from  the 
top  30  feet  of  the  Grassy  Creek  Shale  in  the  deeper  part 
of  the  Illinois  Basin;  however,  the  top  of  the  Grassy 
Creek  was  not  sampled  from  enough  locations  that  it 
can  be  unequivocally  stated  that  the  layer  with  higher 
than  usual  trace  metals  concentrations  is  continuous 
(Frost  et  al.  1985,  p.  69). 

Results  from  the  three  cores  (12IL,  13IL,  and  01IN) 
indicate  that  the  highest  C/S  ratios  and  the  greatest 
enrichments  in  the  metals  are  at  the  top  of  the  Grassy 
Creek.  With  depth,  the  C/S  ratios  decrease  and  metal 
enrichments  decline  such  that  average  values  occur 
several  feet  down  in  the  formation.  Lithologic  features 
(thinly  laminated  beds  and  little  or  no  bioturbation) 
indicate  the  metal-rich  core  samples  are  from  a  com- 
pletely anoxic  environment. 

Nineteen  of  the  22  samples  plot  in  the  high  C/S 
field  (fig.  4B),  and  the  low  level  of  pyrite  is  very  notice- 
able in  the  samples  with  highest  organic  carbon  con- 
tents. The  samples  have  slightly  lower  iron  contents 
than  the  average  for  the  Grassy  Creek  Shale  (table  3), 
and  the  iron  contents  are  independent  of  the  organic 
carbon  contents  (fig.  10A).  The  DOPj  values  indicate  an 
appreciable  amount  of  nonpyritic  iron  in  most  of  the 
samples  (fig.  10B).  However,  DOPr  values  determined 
for  18  samples  show  that  little  of  this  iron  is  reactive 
(fig.  10C);  therefore,  pyrite  formation  was  probably 
iron-limited.  The  anomalously  metal-rich  sample 


(S00127)  of  Hannibal  Shale  from  core  07IL  (noted  in 
table  2  and  included  in  table  3)  has  the  same  charac- 
teristics as  the  Grassy  Creek  samples,  except  that  it  has 
abundant  reactive  iron  (fig.  10C). 

Compared  with  average  Grassy  Creek  Shale,  the  22 
samples  contain  twice  as  much  vanadium  and  certain 
of  the  chalcophilic  trace  elements  (Zn,  Se,  Sb,  and  Ag) 
and  are  slightly  enriched  in  nickel  and  chromium.  They 
generally  contain  the  same  amount  of  the  chalcophilic 
elements  Cu,  As,  Sn,  and  Pb  as  the  average  Grassy  Creek 
Shale  sample  and  lower  amounts  of  cobalt  and  molyb- 
denum than  the  average  Grassy  Creek  Shale  sample. 
Arsenic  and  cobalt  are  commonly  concentrated  in  py- 
rite, so  their  low  concentrations  may  be  controlled  by 
the  low  pyrite  contents.  Because  the  accumulated  ele- 
ments are  chalcophilic,  sphalerite  (ZnS)  and  other  trace 
metal  sulfides  (and  perhaps  antimonides  and  sele- 
nides)  are  likely  to  be  present.  Ripley  et  al.  (1990)  iden- 
tified sphalerite,  chalcopyrite,  nickel-rich  mackinawite, 
and  clausthalite  from  electron  microprobe  analysis  of 
shales  from  the  New  Albany  in  south-central  and 
southeastern  Indiana. 

Lechler  and  Leininger  (1978)  also  reported  for  the 
01  IN  (Sullivan  County)  core  higher  concentrations  of 
zinc,  vanadium,  and  other  trace  metals  in  samples  from 
the  top  of  the  Clegg  Creek  Member  (equivalent  to  the 
upper  part  of  the  Grassy  Creek  in  Indiana).  Higher 
concentrations  of  zinc  and  the  heavy  metals  V,  Ag,  Cd, 
Cr,  Cu,  Mo,  Ni,  and  Pb  were  also  found  at  the  top  of  the 
Clegg  Creek  Member  in  cores  from  Marion,  Brown, 
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and  Clark  Counties  in  south-central  Indiana  (Lechler  et 
al.  1979).  The  grtatest  enrichment  in  trace  metals  found 
in  shale  of  the  New  Albany  in  the  Illinois  Basin,  how- 
ever, was  in  these  cores  in  the  beds  overlying  the  Clegg 
Creek.  For  example,  in  the  top  1  to  2  inches  of  the 
organic  carbon-rich  Henryville  Bed  and  in  the  thin, 
phosphatic  Falling  Run  Bed  (Lechler  et  al.  1979,  Shaffer 
et  al.  1981),  vanadium  concentrations  were  as  high  as 
5,200  ppm  and  zinc  concentrations  were  as  high  as 
17,000  ppm.  Ripley  et  al.  (1990)  also  reported  similarly 
high  metal  concentrations  in  the  Henryville  Bed  in  two 
cores  from  Bartholomew  County,  Indiana. 

Hailer  et  al.  (1982)  described  high  ratios  of  carbon 
to  sulfur,  similar  to  values  reported  above,  in  the  top 
10  ft  of  the  New  Albany  in  cores  from  Sullivan,  Marion, 
Brown,  Jackson,  and  Clark  Counties,  Indiana  (shown  in 
fig.  2  and  plotted  in  fig.  4B).  Hasenmueller  and  Leinin- 
ger  (1987)  state  that  the  low-sulfur  interval  varies  in 
thickness  from  15.5  ft  in  northern  Jackson  County  to  2.5 
ft  in  Clark  County. 

Abundant  evidence  indicates  that  near  the  end  of 
the  deposition  of  the  Grassy  Creek  Shale,  the  rate  of 
sedimentation  was  very  slow  in  the  deep  basin.  Potter 
et  al.  (1980)  found  that  in  slowly  deposited  sediments 
there  is  a  greater  accumulation  of  vanadium  relative  to 
nickel,  presumably  because  organic  matter  long  in  con- 
tact with  seawater  absorbs  more  vanadium.  Hence,  the 
fact  that  the  22  deep  basin  samples  are  much  more 
enriched  in  vanadium  than  in  nickel  suggests  slow 
sedimentation.  Hailer  et  al.  (1982)  related  the  low  ^S 
content  of  the  pyriric  sulfur  in  the  Henryville  and  other 
beds  at  the  top  of  the  New  Albany,  in  Brown,  Clark, 
Harrison,  Jackson,  Marion,  and  Sullivan  Counties,  to  a 
slow  rate  of  clastic  sedimentation  at  the  close  of  the 
deposition  of  the  New  Albany.  From  lithologic  charac- 
teristics of  marker  beds  in  several  cores  and  strati- 
graphic  relations  between  them,  Lineback  (1969) 
concluded  that  by  the  end  of  the  Devonian  or  the  begin- 
ning of  the  Mississippian,  sediment-starved  basin  con- 
ditions prevailed  in  the  deeper  part  of  the  Illinois  Basin. 

Sedimentation  probably  slowed  as  the  sea  in  the 
Illinois  Basin  deepened  (due  to  basin  subsidence)  at  the 
end  of  the  deposition  of  the  Grassy  Creek  Shale.  Be- 
cause the  water  column  was  stratified  and  there  was 
little  circulation  or  oxygenation  of  the  colder,  denser 
bottom  waters,  organic  matter-rich  shales  were  pre- 
served in  the  anaerobic  conditions  at  the  bottom  of  the 
sea.  The  sea  probably  continued  to  deepen  as  the  Han- 
nibal Shale  was  being  deposited  in  southeastern  Illi- 
nois. Eastward  in  Indiana,  little  sediment  was  de- 
posited. Ettensohn  and  Barron  (1982)  suggest  that  the 


sea  depth  eventually  increased  to  the  point  that  up- 
welling  oceanic  waters  off  the  continental  margin  to  the 
south  spilled  over  into  the  epicontinental  Devonian 
Sea;  greater  amounts  of  phosphate  and  heavy  metals 
(in  the  Falling  Run  and  Henryville  Beds),  as  well  as 
much  organic  matter,  accumulated  in  sediments  set- 
tling from  these  nutrient-rich  waters.  Beier  and  Hayes 
(1989)  consider  a  bloom  in  the  production  of  organic 
matter  as  a  principal  cause  for  the  phenomenom  of  the 
Falling  Run  and  Henryville  Beds. 

Blocher  Shale  Formation 

The  relationships  between  carbon,  sulfur,  and  iron  in 
the  New  Albany  Group  formations  discussed  in  pre- 
vious sections  do  not  apply  to  the  calcareous  and 
dolomitic  Blocher  Shale.  Twenty-four  samples  of  Blocher 
Shale  were  analyzed;  they  had  CO2  contents  equivalent 
to  7.4%  to  46.8%  dolomite  plus  calcite.  These  calcula- 
tions assume  that  all  calcium  is  present  in  carbonate 
and  that  the  remainder  of  the  carbonate  carbon  is  bal- 
anced by  magnesium. 

Results  summarized  for  the  Blocher  Shale,  both  on 
the  whole  sample  and  carbonate-free  sample  basis 
(table  2),  show  that  sulfur  content  bears  no  relationship 
to  organic  carbon  content.  Low  iron  and  relatively  high 
and  constant  DOPj  values  indicate  iron-limited  pyrite 
formation. 

Iron  correlates  with  the  detrital  content,  as  esti- 
mated by  clay  and  quartz,  in  the  set  of  samples  (the 
correlation  coefficient  between  iron  and  aluminum  is 
0.62,  significant  at  the  95%  confidence  level).  The  iron 
limitation  arises,  as  expected,  because  iron-rich  clays 
are  diluted  with  iron-poor  carbonate.  For  the  five  Blocher 
samples  analyzed  for  HCl-soluble  iron,  high  DOPr  val- 
ues of  0.73  to  0.96  (fig.  8)  suggest  iron-limited  pyrite 
formation. 

Pyrite  would  be  expected  to  be  iron-limited  in 
nearly  pure  carbonate  samples.  Gibson  (1985),  who 
studied  organic-rich  calcareous  siltstones  of  the  Cam- 
brian Currant  Bush  Limestone  in  Australia,  reported 
that  iron  was  a  major  limiting  factor  in  pyrite  formation 
in  rocks  with  carbonate  as  low  as  32%. 

Organic  carbon  concentrations  in  the  24  samples  of 
Blocher  Shale  analyzed  ranged  from  3.3%  to  9.0%  (5.2% 
to  12.1%  on  the  carbonate-free  basis).  The  consistent 
organic  richness  of  the  Blocher  Shale  and  the  high  DOP 
values  support  the  interpretation  offered  by  Cluff 
(1980)  that  the  Blocher  was  deposited  in  an  anoxic 
deep-water  environment;  but  the  carbon,  sulfur,  and 
iron  relationships  do  not  further  define  the  deposi- 
tional  environment. 


SUMMARY  AND  CONCLUSIONS 


In  374  samples  from  the  New  Albany  Group  in  the 
Illinois  Basin,  mainly  from  within  Illinois,  the  concen- 
tration of  total  sulfur,  which  is  taken  as  a  reliable  meas- 
ure of  pyritic  sulfur,  generally  correlates  with  organic 
carbon  content.  The  relationship  is  consistent  with  the 


theory  that  bacterial  reduction  of  sulfate  concomitant 
with  bacterial  oxidation  of  organic  matter  in  the  sedi- 
ments controls  sulfide  production  and  pyrite  forma- 
tion. An  increase  in  the  production  and  deposition  of 
organic  matter  leads,  particularly  in  deep  water  with 
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restricted  circulation,  to  increasingly  anoxic  conditions. 
Such  conditions  favor  the  growth  and  functioning  of 
anaerobic,  sulfa te-reducing  bacteria,  which,  in  turn, 
govern  the  formation  of  sulfide. 

Conditions  for  preservation  of  organic  material  in 
shale  varied  widely  in  the  Devonian-Mississippian  sea 
in  the  Illinois  Basin.  Lithofacies  data  from  Harvey  et  al. 
(1977),  Cluff  (1980),  and  Cluff  et  al.  (1980)  confirm  that 
progressively  greater  amounts  of  organic  carbon  pre- 
served in  the  shale  correlate  with  progressively  anaero- 
bic conditions.  Such  conditions  were  probably  con- 
trolled by  paleogeography,  and  correlate  with  the  re- 
gional setting  of  the  New  Albany  Shales  proposed  by 
Cluff  (1980). 

The  greenish  gray  Saverton  and  Hannibal  Shales 
were  sampled  mainly  from  the  northwest  part  of  the 
Illinois  Basin.  These  samples  have  organic  carbon 
contents  generally  ranging  from  0.2%  to  2%  and  or- 
ganic carbon  to  sulfur  ratios  averaging  about  1.8  and 
reflect  a  depositional  environment  that  was  generally 
shallow  and  primarily  dysaerobic.  The  degree  of  pyri- 
tization  of  iron  (total  or  reactive)  increases  rapidly  with 
organic  carbon  content;  this  increase  is  an  expected 
consequence  of  the  high  ratio  of  sulfur  to  carbon. 

The  Sweetland  Creek  Shale  and  part  of  the  Grassy 
Creek  Shale  in  western  Illinois  also  are  relatively  low 
(1%  to  3%)  in  organic  carbon  content.  Organic  carbon 
content  in  those  units  generally  increases  southeast- 
ward across  the  shelf  margin,  basin  slope,  and  into  the 
deep  basin  (Frost  1980,  fig.  X-2).  The  brownish  black 
and  black  shales  of  the  Grassy  Creek  and  Blocher  For- 
mations in  southeastern  Illinois  have  organic  carbon 
contents  generally  of  4%  to  9%  and  as  high  as  15%. 
Increasing  carbon  content  indicates  increasingly  anoxic 
depositional  environments.  In  these  samples,  sulfur 
correlates  with  organic  carbon  content,  and  the  average 
ratio  of  organic  carbon  to  sulfur  is  2.8.  DOPT  and  DOPR 
also  generally  correlate  with  organic  carbon  content. 
Variability  in  the  iron  content  of  the  samples  produces 
a  great  deal  of  scatter  in  the  plots  of  DOP  against  organ- 
ic carbon,  especially  as  organic  carbon  content  increases 
to  a  level  at  which  average  iron  content  becomes  the 
limiting  factor  to  pyrite  formation. 

Pyrite  formation  was  found  to  be  iron-limited  in 
the  Blocher  Shale,  which  contains  abundant  carbonate. 


Iron  levels  in  the  shales  were  generally  related  to  the 
clastic  or  detrital  content.  A  small  increase  in  total  iron 
content  as  carbon  content  increases  is  probably  due  to 
adsorption  or  complexation  between  iron  compounds 
and  organic  material.  A  small  increase  in  the  reactive 
iron  content  as  organic  carbon  content  increases  provides 
sufficient  iron  for  pyrite  content  to  increase  and  corre- 
late with  organic  carbon  content  that  is  as  high  as  15%. 

This  study,  because  it  was  able  to  draw  on  informa- 
tion for  a  very  large  number  of  samples,  contributes  to 
a  comphrehensive  picture  of  the  paleoenvironments 
that  existed  during  deposition  of  the  New  Albany 
Group  in  the  Illinois  part  of  the  Illinois  Basin.  But  the 
picture  is  not  complete,  particularly  with  respect  to 
explaining  the  changes  in  conditions  that  led  to 
changes  in  formations  in  different  parts  of  the  Illinois 
Basin. 

No  evidence  from  this  study  substantiates  the  exist- 
ence of  widespread  euxinic  conditions  in  the  Illinois 
Basin  within  Illinois  during  the  deposition  of  the  New 
Albany  Group.  Completely  anoxic  conditions  no  doubt 
prevailed  in  the  deep  basin  in  southeastern  Illinois 
during  deposition  of  finely  laminated  shales,  particu- 
larly of  the  Blocher  and  Grassy  Creek  Shales.  Under 
these  conditions,  some  hydrogen  sulfide  probably  dif- 
fused out  of  those  sediments  into  the  seawater. 

In  contrast,  Beier  and  Hayes  (1989)  and  Ripley  et  al. 
(1990)  concluded  that  formation  of  iron  sulfide  in  a 
hydrogen  sulfide-bearing  water  column  was  signifi- 
cant during  the  deposition  of  the  upper  part  of  the 
Clegg  Creek  Member  and  in  the  overlying  Falling  Run 
and  Henryville  Beds  in  southern  Indiana,  where  depo- 
sition was  on  a  deep  shelf  (fig.  2).  The  same  deep  shelf 
conditions  appear  consistent  with  the  data  (table  2)  for 
the  Bullitt  County,  Kentucky,  core  (02KY).  Interpreta- 
tion of  the  depositional  environment  for  this  core  was 
not  addressed  in  this  report  because  it  is  regionally 
isolated  with  respect  to  the  rest  of  the  data  points.  More 
data  are  needed  from  southwestern  Indiana  and  north- 
western Kentucky  to  document  differences  in  conditions 
between  the  deep  basins  in  those  areas  and  in  south- 
eastern Illinois  and  the  eastern  and  southern  margins  of 
the  Basin  during  Devonian  and  Mississippian  time. 
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Appendix  1  Concentrations  of  total  iron  and  HCI-soluble  iron  in  71  samples 
from  the  New  Albany  Group  determined  by  atomic  absorption  spectro- 
photometry, and  DOPr  calculated  from  pyritic  iron  and  HCI-soluble  iron  concen- 
trations. Samples  are  described  in  Frost  et  al.  (1985). 


Total  Fe* 

HCI-soluble 

Geol.  no. 

Sample  no. 

as  FeO  (%) 

Fe  as  FeO  (%) 

DOPr 

06IL05C1 

S00074 

4.78 

2.34 

0.079 

06IL07C1 

S00076 

3.88 

2.01 

0.18 

06IL09C1 

S00078 

4.81 

2.46 

0.12 

06IL11C1 

S00080 

5.27 

2.09 

0.27 

06IL13C1 

S00082 

4.93 

1.90 

0.27 

06IL15C1 

S00084 

5.10 

1.77 

0.36 

06IL17C1 

S00086 

5.00 

1.48 

0.56 

06IL19C1 

S00088 

5.45 

1.17 

0.66 

06IL21C1 

S00090 

5.64 

1.40 

0.64 

06IL22C1 

S00091 

4.78 

0.944 

0.65 

06IL23C1 

S00092 

5.28 

1.73 

0.55 

06IL24C1 

S00093 

3.39 

1.22 

0.50 

06IL25C1 

S00094 

4.04 

0.815 

0.71 

07IL01L2 

S00127 

4.07 

1.03 

0.32 

11IL02C1 

S00330 

3.71 

0.442 

0.88 

11IL04C1 

S00332 

4.08 

2.13 

0.52 

11IL14C1 

S00342 

4.36 

2.40 

0.46 

11IL15C1 

S00343 

3.78 

1.53 

0.57 

11IL16C1 

S00344 

3.44 

1.82 

0.49 

12IL01C1 

S00402 

5.78 

0.706 

0.83 

12IL01C2 

S00403 

4.28 

0.814 

0.68 

12IL01C3 

S00404 

3.35 

0.525 

0.69 

12IL01C4 

S00405 

4.41 

0.600 

0.79 

13IL04C1 

S00407 

3.60 

I.75 

0.34 

13IL05C1 

S00408 

3.79 

0.652 

0.64 

13IL07C2 

S00410 

3.46 

0.902 

0.67 

13IL08C1 

S00411 

2.72 

0.302 

0.83 

13IL10C1 

S00413 

5.61 

0.626 

0.87 

13IL12C3 

S00415 

6.45 

0.582 

0.89 

13IL14C3 

S00417 

5.24 

1.02 

0.75 

13IL16C2 

S00419 

5.17 

0.988 

0.74 

13IL18C2 

S00421 

5.14 

1.36 

0.65 

13IL19C1 

S00422 

4.42 

1.82 

0.46 

13IL21C2 

S00424 

3.92 

1.02 

0.61 

13IL23C4 

S00426 

4.68 

0.721 

0.80 

13IL26C3 

S00429 

1.74 

0.058 

0.96 

13IL27C1 

S00430 

2.23 

0.077 

0.96 

01IN01C1 

S00352 

3.64 

0.728 

0.67 

01IN01C2 

S00353 

3.50 

0.760 

0.70 

01IN02C1 

S00354 

4.10 

0.836 

0.71 

01IN11C1 

S00372 

5.53 

1.22 

0.74 

01KY07C1 

S00007 

6.03 

1.11 

0.81 

01KY08L1 

S00033 

4.59 

1.96 

0.58 

NAS-0061 

S00152 

4.86 

1.22 

0.61 

NAS-0065 

S00154 

3.88 

0.739 

0.71 

NAS-0093 

S00180 

3.38 

0.511 

0.80 

NAS-0095 

S00182 

4.69 

0.995 

0.69 

NAS-0114 

S00191 

4.48 

1.05 

0.78 

NAS-0119 

S00195 

4.33 

0.604 

0.84 

NAS-0123 

S00198 

4.67 

1.28 

0.62 

NAS-0126 

S00201 

2.80 

0.772 

0.73 

NAS-0129 

S00202 

3.50 

0.866 

0.68 

NAS-0131 

S00204 

3.39 

0.464 

0.85 

NAS-0159 

S00215 

4.04 

0.593 

0.72 

NAS-0166 

S00219 

4.17 

0.779 

0.68 

NAS-0185 

S00227 

5.18 

0.955 

0.70 

NAS-0189 

S00228 

5.31 

0.773 

0.82 

NAS-0208 

S00237 

4.85 

1.68 

0.58 

NAS-0213 

S00239 

3.31 

0.552 

0.75 

NAS-0247 

S00264 

4.99 

0.579 

0.85 

NAS-0249 

S00255 

3.36 

0.484 

0.84 

NAS-0253 

S00257 

3.79 

0.716 

0.74 

NAS-0320 

S00277 

5.58 

1.61 

0.66 

NAS-0343 

S00283 

4.14 

0.771 

0.76 

NAS-0360 

S00289 

5.09 

1.80 

0.58 

NAS-0396 

S00299 

4.06 

1.16 

0.64 

NAS-0421 

S00306 

4.72 

0.949 

0.77 

NAS-0427 

S00308 

3.04 

0.524 

0.77 

NAS-0431 

S00309 

6.06 

3.08 

0.42 

NAS-0447 

S00313 

4.23 

1.01 

0.67 

NAS-0451 

S00314 

4.23 

1.50 

0.55 

'Procedure  Replicate  0.1  g  samples  of  dried  500°C  ash  (used  because 
organic  matter  in  the  whole  shale  interferes  in  a  wet  acid  digestion)  were 
digested  in  hydrofluoric  acid  and  aqua  regia.  Excess  hydrofluoric  acid  was 
complexed  with  boric  acid,  and  the  solution  was  taken  to  50  mL  volume  with 
water  for  measurement  of  iron.  Average  mean  error  between  replicates  was 
0.7%  relative. 
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